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Abstract 



Our representation for solution has the following form 



We present applications of variational - wavelet approach 
to different forms of nonlinear (rational) rms envelope 
equations. We have the representation for beam bunch os- 
cillations as a multiresolution (multiscales) expansion in 
the base of compactly supported wavelet bases. 



1 INTRODUCTION 

In this paper we consider the applications of a new nume- 
rical-analytical technique which is based on the methods 
of local nonlinear harmonic analysis or wavelet analysis to 
the nonlinear root-mean-square (rms) envelope dynamics 
[1]. Such approach may be useful in all models in which it 
is possible and reasonable to reduce all complicated prob- 
lems related with statistical distributions to the problems 
described by systems of nonlinear ordinary/partial differ- 
ential equations. In this paper we consider an approach 
based on the second moments of the distribution functions 
for the calculation of evolution of rms envelope of a beam. 
The rms envelope equations are the most useful for analy- 
sis of the beam self-forces (space-charge) effects and also 
allow to consider both transverse and longitudinal dynam- 
ics of space-charge-dominated relativistic high-brightness 
axisymmetric/asymmetric beams, which under short laser 
pulse-driven radio-frequency photoinjectors have fast tran- 
sition from nonrelativistic to relativistic regime [1]. Anal- 
ysis of halo growth in beams, appeared as result of bunch 
oscillations in the particle-core model, also are based on 
three-dimensional envelope equations [2]. From the for- 
mal point of view we may consider rms envelope equations 
after straightforward transformations to standard Cauchy 
form as a system of nonlinear differential equations which 
are not more than rational (in dynamical variables). Be- 
cause of rational type of nonlinearities we need to con- 
sider some extension of our results from [3]-[10], which are 
based on application of wavelet analysis technique to vari- 
ational formulation of initial nonlinear problems. Wavelet 
analysis is a relatively novel set of mathematical methods, 
which gives us a possibility to work with well-localized 
bases in functional spaces and give for the general type of 
operators (differential, integral, pseudodifferential) in such 
bases the maximum sparse forms. Our approach in this 
paper is based on the generalization [11] of variational- 
wavelet approach from [3]-[10], which allows us to con- 
sider not only polynomial but rational type of nonlineari- 
ties. 
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which corresponds to the full multiresolution expansion in 
all time scales. Formula (|l]) gives us expansion into a slow 
part z^"™ and fast oscillating parts for arbitrary N. So, we 
may move from coarse scales of resolution to the finest one 
for obtaining more detailed information about our dynami- 
cal process. The first term in the RHS of equation (1) corre- 
sponds on the global level of function space decomposition 
to resolution space and the second one to detail space. In 
this way we give contribution to our full solution from each 
scale of resolution or each time scale. The same is correct 
for the contribution to power spectral density (energy spec- 
trum): we can take into account contributions from each 
level/scale of resolution. In part 2 we describe the different 
forms of rms equations. In part 3 we present explicit ana- 
lytical construction for solutions of rms equations from part 
2, which are based on our variational formulation of ini- 
tial dynamical problems and on multiresolution represen- 
tation [11]. We give explicit representation for all dynami- 
cal variables in the base of compactly supported wavelets. 
Our solutions are parametrized by solutions of a number 
of reduced algebraical problems from which one is nonlin- 
ear with the same degree of nonlinearity and the rest are 
the linear problems which correspond to particular method 
of calculation of scalar products of functions from wavelet 
bases and thek derivatives. 

2 RMS EQUATIONS 

Below we consider a number of different forms of RMS en- 
velope equations, which are from the formal point of view 
not more than nonlinear differential equations with ratio- 
nal nonlinearities and variable coefficients. Let f{xi, X2) 
be the distribution function which gives full information 
about noninteracting ensemble of beam particles regard- 
ing to trace space or transverse phase coordinates (a;i, X2). 
Then we may extract the first nontrivial bit of 'dynamical 
information' from the second moments 



x^/(xi,a;2)dxida;2 
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= <xl>=J J xlf{xi,X2)dxidx2 (2) 

Cr^^j,^ = <XiX2>^ J j XiX2f{xi,X2)dxidX2 

RMS emittance ellipse is given by ^ms =< ^1 

x\ > — < X1X2 Expressions for twiss parameters 

are also based on the second moments. 



We will consider the following particular cases of rms 
envelope equations, which described evolution of the mo- 
ments (1) ([1],[2] for full designation): for asymmetric 
beams we have the system of two envelope equations of 



the second order for and a^^ : 



'^I'.+^l.^ + f^I, (^) <^..= (3) 
)7')+eLi/<7', 



The envelope equation for an axisymmetric beam is a par- 
ticular case of preceding equations. 

Also we have related Lawson's equation for evolution 
of the rms envelope in the paraxial limit, which governs 
evolution of cylindrical symmetric envelope under external 
linear focusing channel of strenghts Kr'. 



7' 



Krcr = 



af3'^-f^ ' (73/3272' 



(4) 



where Kr = —Fr/rf3'^jmc^, [i = viy/c ~ •\/l — 7~2 
According [2] we have the following form for envelope 
equations in the model of halo formation by bunch oscil- 
lations: 
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where X(s), Y(s), Z(s) are bunch envelopes, ^x,£,y, = 
F{X,Y,Z). 

After transformations to Cauchy form we can see that all 
this equations from the formal point of view are not more 
than ordinary differential equations with rational nonlin- 
earities and variable coefficients (also,b we may consider 
regimes in which 7, 7' are not fixed functions/constants but 
satisfy some additional differential constraint/equations, 
but this case does not change our general approach). 

3 RATIONAL DYNAMICS 

Our problems may be formulated as the systems of ordi- 
nary differential equations 

Qi{x)'^ ^ Pi{x,t), X = (xi, ...,x„), (6) 
z = l,...,n, max deg Pi — p, maxdegQi — q 

i i 

with fixed initial conditions 2:^(0), where Pi,Qi are not 
more than polynomial functions of dynamical variables xj 



and have arbitrary dependence of time. Because of time di- 
lation we can consider only next time interval: < i < 1. 
Let us consider a set of functions 



^,{t) = x, — {Q,y,) + Piy, 
at 



and a set of functionals 



F,{x) = / <^>^it)dt - Q,x^y^ \l, 
Jo 



(7) 
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where yi{t) {yi{0) — 0) are dual (variational) variables. It 
is obvious that the initial system and the system 



F,{x) = 



(9) 



are equivalent. Of course, we consider such Qi (x) which 
do not lead to the singular problem with Qi (x), when t = 
or t = 1, i.e. Q^ix{0)), Q^{x{l)) ^ 00. 

Now we consider formal expansions for Xi^yf. 

x,{t)^x,{0) + Y,^'iMt) yj{t)^Y.'^Wr{t), (10) 



where (pk{t) are useful basis functions of some functional 
space {L^,LP, Sobolev, etc) corresponding to concrete 
problem and because of initial conditions we need only 

^fe(0) = 0,r = l,...,7V, i = l,...,n, 



A = {Aa = {Aa = (Ai,A2,...,Af) 



(11) 



where the lower index i corresponds to expansion of dy- 
namical variable with index i, i.e. Xi and the upper index 
r corresponds to the numbers of terms in the expansion of 
dynamical variables in the formal series. Then we put ( |l0| ) 
into the functional equations (H) and as result we have the 
following reduced algebraical system of equations on the 
set of unknown coefficients A^' of expansions (|0|): 



L{Q,j,X,ai) = M{P,j,X,pj), 



(12) 



where operators L and M are algebraization of RHS and 
LHS of initial problem where A (|rT]) are unknowns of 
reduced system of algebraical equations (RSAE)(p^. 

Qij are coefficients (with possible time dependence) of 
LHS of initial system of differential equations (||) and as 
consequence are coefficients of RSAE. 

Pij are coefficients (with possible time dependence) of 
RHS of initial system of differential equations (||) and as 
consequence are coefficients of RSAE. / = (ii, ...,1^+2), 
J = {ji, jp+i) are multiindexes, by which are labelled 
aj and /3/ — other coefficients of RSAE (113): 



i<i/c<P+i 

where p is the degree of polinomial operator P i 



(14) 



where q is the degree of polynomial operator Q (g), = 
(1, + 2), ip'i^ = dcpijdt. 

Now, when we solve RS AE ( [T^ and determine unknown 
coefficients from formal expansion ( p^ we therefore ob- 
tain the solution of our initial problem. It should be noted 
if we consider only truncated expansion ( p^ ) with N terms 
then we have from ( p^ the system of x n algebraical 
equations with degree £ = max{p, q} and the degree of 
this algebraical system coincides with degree of initial dif- 
ferential system. So, we have the solution of the initial 
nonlinear (rational) problem in the form 



^'•j\/'/\/^,/-j'Af'/^^A''ipMf'M^V-^^ 
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where coefficients \\ are roots of the corresponding re- 
duced algebraical (polynomial) problem RSAE ([1^. Con- 
sequently, we have a parametrization of solution of initial 
problem by solution of reduced algebraical problem (|l2|). 
The first main problem is a problem of computations of 
coefficients a/ (p^, (3j (13) of reduced algebraical sys- 
tem. These problems may be explicitly solved in wavelet 
approach. The obtained solutions are given in the form 
(p3|), where Xk{t) are basis functions and are roots 
of reduced system of equations. In our case Xk{t) are 
obtained via multiresolution expansions and represented 
by compactly supported wavelets and \\ are the roots of 
corresponding general polynomial system (|l2|). Our con- 
structions are based on multiresolution approach. Because 
affine group of translation and dilations is inside the ap- 
proach, this method resembles the action of a microscope. 
We have contribution to final result from each scale of res- 
olution from the whole infinite scale of spaces. More ex- 
actly, the closed subspace Vj{i G Z) corresponds to level 
j of resolution, or to scale j. We consider a multiresolution 
analysis of L^(R") (of course, we may consider any dif- 
ferent functional space) which is a sequence of increasing 
closed subspaces Vf ...V-2 C V-i C Vb C C V"2 C ... 
satisfying the following properties: 

So, on Fig.l we present contributions to bunch oscillations 
from first 5 scales or levels of resolution. It should be noted 
that such representations (1), (15) for solutions of equa- 
tions (3)-(5) give the best possible localization properties 
in corresponding phase space. This is especially important 
because our dynamical variables corresponds to moments 
of ensemble of beam particles. 

In contrast with different approaches formulae (1), (15) 
do not use perturbation technique or linearization proce- 
dures and represent bunch oscillations via generalized non- 
linear localized eigenmodes expansion. 
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Figure 1 : Contributions to bunch oscillations: from scale 

2^ to 2^ 
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